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Abstract. - We study theoretically a scheme in which particles from an incident beam are 
trapped in a potential well when colliding with particles already present in the well. The balance 
between the arrival of new particles and the evaporation of particles from the trapped cloud 
leads to a steady-state that we characterize in terms of particle number and temperature. For 
a cigar shaped potential, different longitudinal and transverse evaporation thresholds can be 
chosen. We show that a resonance occur when the transverse evaporation threshold coincides 
with the energy of the incident particles. It leads to a dramatic increase in phase space density 
with respect to the incident beam. 



The trapping of atomic particles has been an invaluable tool for recent developments 
in atomic physics and quantum optics. It can be performed by suddenly switching on a 
confining potential when the particles are in the vicinity of its minimum. This method is used 
successfully to trap ions in electromagnetic traps Q| , and neutrons Q| or atoms Q in magneto- 
static traps. Another way of trapping particles is to take advantage of a dissipative mechanism, 
such as in a magneto-optical trap Q . In this way, the trap can be loaded continuously since 
a friction force dissipates the excess energy of the particles and prevents them from escaping. 
The loading of neutral atoms into a trap using collisions with a buffer gas also belongs to the 
dissipative category 

In this paper we investigate a different loading mechanism, in which particles are injected 
into a potential well, and can be trapped by undergoing an elastic collision with one of the 
particles already present in the well. After this collision, the incident particle has an energy 
below the depth U of the potential well and gets trapped. The excess energy is then redis- 
tributed to the whole trapped sample via elastic collisions and it is subsequently released by 
the evaporation of a trapped particle. We show that it is possible to accumulate in this way 
particles with an equilibrium temperature T <C U/k^. We also show that, for experimentally 
feasible conditions, a large increase in phase space density can be achieved with respect to 
the one of the incident beam of particles. Our study constitutes a realistic description of the 
continuous loading of a trap, consistent with Liouville's theorem It has to be contrasted 
with the elegant, but simplified model of [Q, where new atoms are injected into a trap at a 
given energy, but where evaporation of previously trapped atoms at the same energy is not 
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Fig. 1 - (a) First model; atoms with an energy (1 + e)(7 are injected into a harmonic trap and can be 
captured by collisions with particles already present in the trap. Trapped atoms with an energy larger 
than U are evaporated, (b) Second model; incoming atoms undergo collisions with trapped particles 
or with those reflected on the potential barrier at z < 0. Evaporation occurs if the longitudinal or 
transverse energy of a particle exceeds the evaporation threshold Uz or U±. The phase space density of 
the trapped particles can be orders of magnitude larger than the phase space density of the incoming 
beam. 



taken into account; in that case, the gain in phase space is only hmited by ad-hoc losses, such 
as three-body colHsions. 

We develop two models in this paper. The first model deals with an isotropic harmonic 
trap, assuming an evaporation criterion based on the total energy of a particle. This allows to 
use simple analytical expressions for the evaporation rate and it leads to a proof-of-principle 
of the process, as long as the energy of the injected particles is small with respect to the trap 
depth. The second model is more elaborate and leads to much more spectacular predictions 
for a realistic situation. It assumes an anisotropic trap (cigar shaped potential) and it takes 
advantage of different evaporation rates along the long axis of the cigar and in the transverse 
directions. One can then accumulate particles in the well even if their incident energy notably 
exceeds the well depth. We conclude the paper by giving some indications on the dynamics 
of the loading process. 

For our first model, we consider an isotropic harmonic potential with frequency lu. An 
atomic beam with flux $ and a mean energy U{1 + e) per particle is injected into a poten- 
tial weU (fig. 0a). We assume that those atoms undergoing collisions are trapped. Atoms 
are evaporated as soon as their total (kinetic-|-potential) energy after a collision exceeds the 
threshold U. The steady-state energy distribution P{E) of the trapped gas is approximated 
by a truncated Boltzmann distribution with temperature T [||: 

P{E) (X e"-^/"^ e{U - E) , 

where 9 is the Heaviside step function. For simplicity we set the Boltzmann constant fce = 1. 

In order to determine T and the number N of trapped atoms, we equalize the incoming 
and outgoing fluxes of energy and particles. An incoming atom brings the energy [/(I -f e). 
The average energy carried away by an evaporated atom can be written U + k T, where n 
is a dimensionless coefficient which depends on the ratio r/ — U/T. This coefficient can be 
obtained by a simple integration over the Boltzmann equation [|[^. From the equation 

{l + e)U = U + kT ^ , (1) 
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we deduce i], hence the temperature. The function K{ri)/ri equals 1/6 for rj — Q and decreases 
when rj increases. Therefore, a steady-state exists only for values of e smaller than 1/6. For 
large 77, the function K(ry) tends to 1 and we obtain ~ e~^, i.e. 

e < 1 ^ T ~ £[/ . (2) 

In other words, the equilibrium temperature is equal to the excess energy of the incoming 
atoms. 

We now determine the number of atoms in steady-state by comparing the flux of atoms 
injected into and evaporated from the trap. We first note that only a fraction of the incident 
fiux <I> contributes to the feeding of the trap. Indeed, the above description of the evaporation 
is only valid in the collision-less (or Knudsen) regime, where the mean free path is much 
larger than the size of the atom cloud. Therefore, most of the incoming atoms cross the trap 
without any collision. To calculate the fraction / of incoming atoms which collide with a 
trapped atom, we assume that all incoming trajectories pass through the center of the trap. 
For the low temperatures considered here, the elastic collisions are isotropic (s-wave regime) 
and they are characterized by the total elastic cross-section a. We find / = \/27rfno(T, 
where f — ^ kT jmuP' is the size of the trapped cloud and no the density at the center of the 
harmonic trap. The outgoing flux can be calculated using classical kinetic equations. The 
average collision rate is 7 = ■\/2/7r ng o Au, where Aw = ^JkT /m. In the collision-less regime 
7 <C w, the probability that an atom is evaporated after a coUision is p ~ 2 776"'' if 77 1 
The balance between the incoming flux $in = /$ and the outgoing fiux <I>out = PlN then 
gives: 

= p-yN ^ N = e- e^/" . (3) 

This result is intuitive. In the regime ~ 77 ^ 1, the number of atoms in the trap 
scales as e^^^ ~ e^/^ , which is a direct consequence of Boltzmann's law. The number N 
can be in principle very large. However one should keep in mind the validity criterion for 
the Knudsen regime / < 1, which imposes e^/'^ < U / {^mauj). This puts an upper bound 
-^max ^ fXJ / {muj^ a) on the number of trapped atoms. 

Qualitatively new features arise in the second model that we now discuss. Here, atoms 
are injected across the plane z — zq with a negative velocity along the z axis, and with some 
transverse potential and kinetic energy (flg. |^b). For such a geometry, the evaporation can 
be either longitudinal or transverse. The longitudinal evaporation is a direct consequence of 
the loading mechanism: an atom which crosses the plane z = zq with a positive velocity is 
evaporated. Furthermore one can experimentally arrange that an atom is also evaporated 
when its distance r = {x^ + y^)^!"^ to the z axis exceeds a threshold value tq. For instance, 
one can use a magnetic gradient to conflne the atoms transversally and use a radio-frequency 
wave to flip their magnetic moment when they cross the surface of the cylinder r — |]ll| . 

The possibility to control independently the transverse and longitudinal evaporation thresh- 
olds is of particular interest if one considers an anisotropic trap with lOz ^ where ujx_ 
and stand for the oscillation frequencies in the xy plane and along the z axis, respectively. 
In this case, one can reach a high-flux regime in which the z motion is in the hydrodynamic 
regime, i.e. 7 ^ All incident atoms undergo a collision with the trapped atoms and are 
captured. We now study how the simple approach presented above is modified and we derive 
the expected gain in phase space density for a realistic situation. 

The evaporation rate along the z axis is notably reduced as compared with the rate derived 
in the collision-less regime for the same ratio 77^ — Uz/T. Indeed most atoms which emerge 
from a collision with an energy larger than the evaporation energy = moj'^ZQ/2 undergo, 
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before reaching the point zq, another coUision which can bring their energy below Uz. Using 
a molecular dynamics simulation jl^-p^, we have calculated (i) the probability Pz that an 
atom reaches z = zq after a collision and is evaporated, (ii) the average energy Uz + Kz T 
carried away by this atom. For rjz in the range from 4 to 7, we find: 

< 7 : Pz- 0.14 e"''^ uJz/j Kz ~ 2.9 . (4) 

Under these operating conditions, the evaporation threshold U± = mu;^rg/2 in the xy 
plane is a crucial control parameter. One can set U± to a value larger than Uz^ so that the 
particles evaporated transversely carry an energy U±_ + k±T notably larger than the energy 
Uz + KzT for a particle evaporated along z. If the trap were operated in the collision-less 
regime for the z motion, this would make the transverse evaporation unlikely and inefficient. 
However, thanks to the reduction of pz by a factor 1^2/7 with respect to the collision-less case, 
we can raise U± to a, value larger than Uz and still obtain similar probabilities pz and p± for 
longitudinal and transverse evaporation. 

The equilibrium temperature T is now obtained by generalizing (|l|): 

il + e)Uz= {Uz + KzT)+ iU^+Ki_T). (5) 

Pz +P± Pz+ P± 

Similarly the number of atoms in steady state is obtained from the generalization of (^: 

^in = ip±+PzhN . (6) 

The values of Pz and Kz have been given in (^. We have also derived the corresponding 
values for the transverse motion from a molecular dynamics simulation. For collision dynamics 
ranging from the Knudsen regime (7 <C to the hydrodynamic regime (up to 7 = 5 lo±), 
we fit our results for ry^ = 8 to 13 by the formula: 

p^ ~ 2.0 e-"^ , K± ~ 2.0 . (7) 

ijj± + 1.4 7 

The set of equations (|4[Q) now allows us to determine T and N for a given experimental 
situation. 

As a concrete example we consider a beam of rubidium atoms (cr — 7.1 x 10^^^ m^) 
injected in a trap with uj± = 100 and Uz/'^t^ = 10 Hz. The flux is $ = 10^ atoms/s, 
with an average velocity u; — 20 cm/s and a velocity dispersion Avi — 4 cni/s along each 
of the three axes x,y,z (temperature Ti ~ 17 ^K) ||l5[| . This beam is confined transversely 
by a harmonic potential with the same frequency un_ as in the trap. In these conditions, the 
initial phase space density is 2?^ ~ 2 x 10^^. The barrier height Uz in zq (see fig. lb) is 
chosen such that only atoms with an incident velocity larger than Vi — Avi are transmitted 
{Uz = m{v, - Avifjl - 135 ^K). The flux passing va. z = z^^ is then <i>in = 0.84 $ and the 
average excess energy tUz of the incoming atoms is ~ 140 /iK (i.e. e ~ 1). 

The temperature and the number of atoms are shown in fig. g as a function of the transverse 
evaporation threshold. We also give the phase space density T) in the trap calculated using 
Boltzmann statistics. It shows a sharp maximum T) ~ 10~^ when U ^ is close to 270 /iK, in 
which case p±/pz ^ 2. The gain in phase space density with respect to the incident beam 
is then ^ 500, while it is only 10 if one restricts evaporation to the longitudinal direction 
(/7j_ = 00). In this optimum regime, we find in steady-state 7 ^ cjj_, N ~ 10^ atoms and 
T = 24 //K [rjz ^ 5.7, yy^ ^ 11.4). Note that the temperature is now much lower than the 
excess energy eUz- This has to be contrasted with the result of the first model where one 
always gets T > eU. 
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Fig. 2 - Variation of the temperature (a), of the number of trapped atoms (b), and of the phase 
space density (c), as a function of the transverse evaporation threshold U±. The characteristics of 
the injected atomic beam are given in the text. The dotted lines in (a) indicate the variation of the 
temperature if one neglects either longitudinal evaporation (line with negative slope) or transverse 
evaporation (horizontal line). 



The spectacular resonance in the steady-state phase space density shown in fig. ^ is a 
general feature of our second model. It occurs when the transverse evaporation threshold 
Uji_ is set to a value near the total incident energy (1 + e)Uz- This resonance, which allows 
to produce quantum degenerate trapped clouds from incident atomic beams still far from 
degeneracy, constitutes the main result of this letter. 

We now discuss the characteristics of the system on the two sides of the resonance. When 
U±^ is lower than the optimum value, the evaporation is essentially transverse since pj_ ^ Pz- 
In this case, (^ simplifies to give 

P± » P. : T = — ((1 + €)Uz - U±)) . (8) 

K± 

We recover in this way the linear variation of T for a small U± (see fig. ||). The number of 
atoms, given by $in ~ p±jN, increases as U± increases. The approximation (|^) is valid until 
U± approaches (1 + e)Uz, in which case it would lead to T ~ 0. At this point the longitudinal 
evaporation cannot be neglected anymore, and Pz and p± are of the same order. 

When U± is larger than the optimal value, the transverse evaporation becomes inefficient 
and one is left with a purely longitudinal evaporation. The temperature and the number of 
trapped atoms are then independent of U±. From (||-^), one gets T = tUz/nz (horizontal 
dashed line in fig. ||a) and N = e''^/(0.14 w^), which is very reminiscent of the results of 
the first model. 

The gain G in phase space is an important result of our study. It depends on the incident 
flux $ for a given loading and trapping configuration (parameters e and w^/wz). For relatively 
small incident flux, we flnd that G increases with increasing We consider again the above 
example and we suppose that we increase $ from 10^ s~^ to 10* s""'^. Due to an increase of 
the evaporation efficiency, we find that the steady-state atom number N is multiplied by 15, 
and that the temperature decreases by 10%. The phase space gain then passes from 500 to 
900. 

The phase space gain G as a function of the fiux $ saturates when the transverse motion 
enters the hydrodynamic regime (7 tj^). In this case, the ratio p^jpz obtained from (|^0) 
is a function of r]}_ — rjz and e only. The temperature and the ratio N/ $ are independent of 
the flux $, hence a gain G also independent of <&. Assuming e ~ 1 as in the example above. 
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this maximal gain is G ~ lO'^ for u)±/ujz = lO'^, and G ^ 10'' for lo^/lUz — 

We now investigate the dynamics of the system and the time necessary to reach the steady- 
state of the system around the resonance. The time evolution of the atom number N and of 
the total energy 3 NT is given by: 

dN 



dt 

d{3NT) 
Jt 



= '^in-iPz+P^hN (9) 

= <Pinil + e)U, - p,^ N {U, + K,T) - p^j N {U^ + K^T) . (10) 



The steady-state of this set of non- linear equations corresponds to the solution of (|[^). Ini- 
tially, the description in terms of a thermal equilibrium fails since no atoms are present. The 
loading of the trap is initiated by collisions between atoms of the incoming and the reflected 
outgoing beams. In the above example, the collision rate 7 between the two beams is of the 
order of ujz, which ensures an efficient start of the whole process. 

The detailed study of the non- linear dynamics involved in (^-10) is outside the scope 



of this letter. A simple hint on the behavior of the solutions consists in linearizing these 
equations around the steady-state discussed above. Suppose for simplicity that evaporation 
occurs essentially in the longitudinal direction. The two time constants of the corresponding 
system are (within a numerical coefficient) ti ~ riql and T2 = t/Vzj where r — Ng/^in is 
the time required to send in the trap the number of atoms corresponding to the steady-state 
number Ng. Since 77^ ^ 1, these two time constants are very different from each other and 
the time ti needed to reach equilibrium exceeds by a factor 77^ the "natural" time scale r. 

For the same experimental conditions as in fig. |^, we find a slow time evolution of N{t), 
with a time constant ri ~ 160 s (to be compared with r = iVs/$in ~ 10 s). By contrast, 
the temperature T{t) reaches rapidly a value close to steady-state. We have checked that a 
numerical simulation based on molecular dynamics gives values for the evolution of N{t) and 
T{t) which differ by no more than 20% from those deduced from the set of equations (PHlOl). 
It is worth noting that for such long time scales losses, e.g. due to 3-body collisions, become 
important and will limit in practice the gain in phase space density as in ||^. 

To summarize, we have shown that a continuous loading of a non-dissipative trap could 
be achieved with the help of evaporation. This process can be remarkably efficient in the case 
of an anisotropic geometry. In this case, the incident particles can have an incident energy 
eUz of the order of the entrance barrier Uz and still reach a temperature T <^ Uz- For a 
realistic configuration a gain in phase space of several orders of magnitude is possible between 
the incident beam and the trapped cloud. This method can therefore be used to produce in 
a continuous way a Bose-Einstein condensate. It constitutes an alternative to the proposal 
of based on a continuous evaporative cooling of a guided beam. It is also complementary 
of the scheme recently achieved at MIT where a succession of condensates was merged in the 
same trap, producing thus a steady-state BEC in a given location n%. 
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